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The complete symmetry characterization of eigenstates in bare opal systems is obtained by means 
of group theory. This symmetry assignment has allowed us to identify several bands that cannot 
couple with an incident external plane wave. Our prediction is supported by layer-KKR calculations, 
which are also performed: the coupling coefficients between bulk modes and externally excited field 
tend to zero when symmetry properties mismatch. 
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I. INTRODUCTION 

Periodic dielectric structures (the so called photonic 
crystals) are probably one of the most exciting topics in 
contemporary physics, not only for its scientific relevance 
but also because of the remarkable applications that have 
been suggested, specially for systems with gaps in the 
visible and infxaredjrange. So, since the first proposal 
in the late eighties ,Era many efforts have been devoted to 
the fabrication of materials with an absolute frequency 
gap in such region. However, material engineering of 3D 
periodic systems at the relevant length scale (hundreds of 
nm to 1 micron) presents non trivial technological prob- 
lems in a straightforward manufacture. Several different 
self-assembly approaches have been proposedH and one of 
the most promising strategies is that based on artificial 
opal systems .ai3 

It is clear that one of the most important aspects of 
sample characterization relies on a proper understanding 
of optical measurements. The structures within the re- 
flectance or transmittance spectra are currently related 
with gaps and pseudogaps in the infinite band structure, 
but the symmetry properties of its eigenfunctions can 
also produce unexpected patterns: symmetry-inactive 
modes (ie. bands that cannot be coupled with incident 
light due to symmetry reasons) have been reported of 
some 2D photonic crystalsaa and some authors have pre- 
dicted the i^ajHe phenomenon to take place in 3D materi- 
als as well.Bli3 Thus, the classification of eigenmodes ac- 
cording to their symmetry properties can provide us valu- 
able information when comparing our calculation with 
experimental results. 

The aim of this work is to present a complete symmetry 
characterization of eigenstates for the abovementioned 
opal structure. In order to reach this goal, the paper is 
organized as follows: In Section II we introduce the solu- 



tion of Maxwell equations as an eigenvalue problem that 
can be restricted to a band structure calculation. Section 
III describes the symmetry properties of electric field ac- 
cording to group theory. In Section IV we present the 
mode assignment for a bare opal system. The existence 
of uncoupled modes is discussed in Section V. Finally, in 
Section VI we summarize our work. 



II. MAXWELL EQUATIONS AND PHOTONIC 
BAND STRUCTURE 

As far as the point at stake is the propagation of light 
in a dielectric material, we shall concern ourselves with 
the macroscopic Maxwell equations. |— If we impose our 
problem to fulfil some requirements,!!-!! we can decouple 
the equations in order to finally arrive to an expression 
entirely in terms of the magnetic field H(r) 



V x ( -4-tV x H(r) 

\e(r) 



-) H(r) 



(1) 



With the definition of a suitable differential operator 



8H(r) = V x ( ——\7 x H(r) 



(2) 



equation ([!]) can be explicitly written like an ordinary 
eigenvalue problem: 



9H(r) 



-) H(r) 



(3) 



The somehow arbitray definition of ensures the oper- 
ator be hermitian, which implies that its eigenfunctions 
have several properties that are extremely useful in nu- 
merical calculations.E-3 
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In the alternate approach, we can eliminate H(r) to 
obtain an equivalent formulation for the stationally elec- 
tric field. This way of considering the problem results 
in a generalized eigenvalue equation with the dielectric 
function playing the role of a new hermitian operator in 
the right side: 

S x E(r) eee V x (V x E(r)) = 
= (^) 2 £ (r)E(r)EE(^) 2 S 2 E(r) (4) 

It is clear that both choices provide the correct physics, 
so the preference for (|J) instead of (|J) (or viceversa) in 
the numerical evaluations will only depend on our per- 
sonal convenience. 

We finally want to remark that Bloch's theorem can be 
applied to photonic crystals due to the periodicity of the 
dielectric function s(r). Thus, we can restrict the eigen- 
value problem to the calculation of the band structure of 
our system and therefore rewrite Equations (|J) and (jij) 
as 

§HJ(r)=^) 2 H^(r) (5) 

g lE J(r) = (^)' 'saEE(r) (6) 

where u> is the frequency of the eigenstates, k the Bloch 
wave vector within the first Brillouin zone and n a dis- 
crete index for increasing eigenfrequencies. 

III. EIGENMODE CHARACTERIZATION 

A. Eigenfunctions and symmetry operators 

Bloch's theorem states that the electric field within a 
photonic crystal can be written as 

E£(r) = e ik ' r u k (r) (7) 

in which k is the Bloch wave vector inside the first Bril- 
louin zone and Uk(r) a periodic vector function of the 
lattice structure. It follows from the abovementioned pe- 
riodicity that 

Ek -e lk r ^f q e tq r (8) 
q 

where q is a reciprocal lattice vector of the structure. 

Let A be a symmetry operator of the lattice. The way 
in which a general symmetry transformation A operates 
over a vector field f(r) is: 

A.f(r) = Af(A _1 r) (9) 



Hence 

A • E£ = Ae <k A " lr f q e iq ' A ~ lr (10) 
q 

Taking the orthogonality of symmetry matrices into ac- 
count, 

A.Eg = Ae lAk r ^ f q e lAq r (11) 
q 

For a wave vector k, the k- vector point group Gk is de- 
fined as the set of symmetry operations which satisfy 

Ak = k + q' (12) 

with q' in the reciprocal lattice. Thus, A being part of 
Gk, 

A • E£ = e i ( k+q '}- r Y Af q e 4Aq r (13) 
q 

Given that the extra phase can always be included in 
the infinite summatory, we conclude that the symmetry 
transformation of E£(r) generates another Bloch func- 
tion with the same wave vector: 

l.E£ = e ikr v k (r) (14) 

Let us now consider the action of a symmetry operator 
A on both sides of Equation (^) : 

(l.H 1 )E^(r)= (^) 2 (A.S 2 )E£(r) (15) 

For any symmetry operator that leaves e(r) invariant 
within the unit cell, it can be easily proved that it com- 
mutes with Si and S 2 . Therefore 

S 1 (l.EL l (r))= (^) 2 § 2 (A.Ej;(r)) (16) 

which implies both E£(r) and (A • E£(r)) be eigenfunc- 
tions with the same eigenfrequency. 

On the assumption that Gk = {Ai} is the point group 
of some wave vector k, it follows from Equations ( |l4| ) and 
(|l|) that E£(r) must verify 

A i E^ g (r) = ^a 1 g!g ,E^ g ,(r) (17) 
g' 

where g is eigenmode degeneracy. We shall henceforth 
concern ourselves with the meaning of that set of scalar 
coefficients {a* g ,}. 
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B. Theory of representations 

Within the framework of the usual group theoryjl! it 
can be straightforwardly proved that the scalar coefficient 
matrices a 1 in ( J17| ) form a representation of (i.e. are 
homomorphic with) Gk group. On the other hand, in 
the case where the degeneracy is said to be normal (i.e. 
when it is not possible to remove such a degeneracy by 
changing a parameter which does not alter the stated 
symmetry of the problem), it may be seen that the set 
{E£ g ( r )} spans an irreducible invariant space under the 
Gk group of symmetry operators. Thus, the {a 1 } set of 
matrixes constitutes an irreducible representation R of 
Gk point group, and R label can be therefore assigned to 
nth eigenmode at k point. 

If Ai belongs to a class G, then 

5>i..=x8 (18) 

9 

with Xc t ne character of class G in some irreducible rep- 
resentation R. Given that each R has its own characters, 
we can use them to distinguish the irreducible represent 
tations in order to classify the different eigenmodes .EJ3 
Further details about character evaluation can be found 
in Appendix A. 



B) with the corresponding modes at adjacent wave vec- 
tors by following the loss of symmetry as we go from one 
point group to another which is a subgroup of the first. 

In Figures 2 to 4 we present the photonic band struc- 
ture of a face-centered cubic lattice of spheres. Results 
are plotted in terms of reduced frequency a/A where a 
and A denote the lattice parameter and the vacuum wave- 
length respectively. The dielectric constants of the sphere 
and the background were assumed to be 2.104 (silica) 
and 1.0 (air). The ratio of the lattice parameter and the 
sphere radius was chosen to verify close-packing condi- 
tion. This band calculation .was carried out by means of 
an iterative implementationE3 of the plane wave expan- 
sion method on Equation (^). Of course, the photonic 
band structure of a bare opal system has been previously 
calculated, c3 so the main outcome of our present work is 
obviously the mode assignment for the electric field. The 
labels within the figures have been determined by the 
numerical evaluation of only several eigenmodes, as far 
as the abovementioned compatibility relations allowed us 
to easily connect the irreducible representations for adja- 
cent wave vectors. For the sake of completeness, we also 
list in Table II the symmetry characterization of eigen- 
modes at r, X, L, U, W and K points for increasing 
values of the eigenfrequcncy. 

V. UNCOUPLED BANDS 



IV. SYMMETRY CHARACTERIZATION IN 

BARE OPAL SYSTEMS A. Symmetry assignment and uncoupled modes 



Synthetic bare opals are constituted by Si02 spheres 
that organize themselves to form a face-centered cubic 
(fee) lattice. As far as this process has been extensively 
described in previous workso, we just mention that opa- 
line sample growth method is based on natural sedimen- 
tation of silica nanospheres along the (111) direction of 
fee lattice. Although the refractive index contrast be- 
tween Si02 and air does not allow bare opals to exhibit 
complete gaps, they can be used a&-,templates to ob- 
tain the so called inverted structures jL3~E2l which are one 
of the most promising strategies between the different 
self-assembly approaches to photonic crystal fabrication. 
Thus, no one can deny that bare opals still constitute a 
topic of present interest. 

In order to work out the classification of eigenstates 
described in the previous section, we shall first concern 
ourselves with the intrinsic properties of the face-centered 
cubic lattice, as far as the symmetry point group which 
determines the set of allowed irreducible representations 
for any wave vector k is irrespective of the particular re- 
alization of our system. Table I lists the symmetry point 
groups at the corners of the irreducible part of the first 
Brillouin zone (see Figure 1) according to Schocnflies no- 
tation. Given that R representation can be assigned to 
nth mode at one of those particular points, we can then 
easily derive all its compatibility relations (see Appendix 



As mentioned in Sec. I, light attenuation produced by 
uncoupled bands ( i.e., bands with eigenmodes that can- 
not be excited by the incident light because of symmetry 
reasons) has been reported of some 2D photonic crys- 
tals. Moreover, some authors have predicted a, similar 
phenomenon to also take place in 3D materials |f~li3 but 
as far as we know, there is no experimental evidence of 
this mechanism operating in an opaline structure. The 
exposition of how these uncoupled modes can be identi- 
fied as a simple consecuence of the symmetry assignment 
will be the aim of the present section. 

In order to reach this goal, let us focus our interest on 
what happens when light incidence is normal to the {111} 
set of planes in the opal. This is the direction along the 
structure naturally grows and from which experimental 
data are usually obtained. If we consider the z axis in 
the normal direction to the sample surface, the incident 
plane waves can be expressed in terms of the following 
vector basis 

{vi = (l,0,0)e lfcz ; v 2 = (0,l,0)e fc } (19) 

where k is the wave number of the plane waves in free 
space. On the other hand, it is clear that vectors in ( JTg| ) 
also form a basis set of some reducible or irreducible rep- 
resentation of C% v point group. In order to determine 
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that representation, we have to apply the symmetry op- 
erations of Czv over the vector basis: 




^v 2 



Ct,vi = -Vi 
5»v 2 = v 2 



(20) 



(21) 



From Eqs. @ and (|l|) we can obtain the trace (char- 
acter) of the transformation matrixes: 



x(C' 3z ) = -l; x(ff») = « 



(22) 



If we compare these results with the properties of C$ v 
point group, we find that these are the characters of the E 
representation. But the boundary condition of tangential 
components of EM field being continous at the interface 
imposes the representation (i.e. the symmetry proper- 
ties) both inside and outside the system to be the same. 
Consequently, only _E-labelled states will contribute to 
light transport along the (111) direction of the structure. 
This kind of restriction in available bands is present in 
several directions of the lattice, but we are mainly inter- 
ested on those which include the T point, i.e. those of 
which some observable consecuences are expected. 

In Table III we summarize the information about un- 
coupled modes along the relevant directions of a face- 
centered cubic structure provided by symmetry consid- 
erations. Let there be a frequency range where only un- 
coupled modes exist. Despite of its nonzero density of 
states, we would expect a total reflection (zero transmis- 
sion) in such a range as a consequence of photonic states 
not being coupled with the incidence plane wave. 

Unfortunately, this is not the case for bare opal sys- 
tems: as shown in Figures 2 to 4, uncoupled modes al- 
ways coexist with available ones, and consequently trans- 
mittance measurements cannot confirm the predictions 
of group theory. Therefore, an alternative test will be 
required. 



B. Estimation of the coupling by means of the 
layer-KKR method 

It is well known that the photonic band structure 
of a regular array of non-overlapping spheres in a host 
medium can be very accurately calculated within the 
framework of the layer-KKR methodJ'EirE£l. Moreover, 
we can also obtain the transmission and reflection coef- 
ficients through a finite slab of such material combining 
the scattering matrices of its different layers. On the 
assumption that slab thickness is very large, the electro- 
magnetic field in the middle of that slab can be regarded 
as a superposition of the true bulk eigenstates. Thus, a 
kind of overlapping integral of this field and the actual 



eigenstates of the infinite system would provide an es- 
timation of the coupling between the structure and the 
incident light. 

Let the z axis be in the normal direction to the slab 
surface. Hence, the electric field between two consecutive 
sphere layers can be expressed in terms of its Fourier 
components as 



E 



void 



(r)=5>+e^- + u-e iK *-- r 



K t = k ll 



q ± z 



"(kll 



(23) 



(24) 



Here, q denotes the 2D reciprocal lattice vector associ- 
ated to the surface periodicity, and fcii = {k x ,k y ) the 
wave vector within the surface Brillouin zone. On the 
other hand, these Fourier components also verify 



U2 



1 = 



(25) 



where Q is an appropriately constructed tranfer matrix 
which relates the electric field in the void between the 
(N — l)-th and JV-th layers to that between the iV-th 
and (N + l)-th ones. Since Q is not hermitian, its left 
eigenstate 



u L = 



(26) 



is generally different from the complex conjugate of the 
right one, but we assume that both are normalized to 
meet the orthogonality condition: 

(u n Jui=6 nn> . (27) 

We can then define the coupling coefficient c„ as 

c n = (O^buik, (28) 

with itbuik a column vector which contains the Fourier co- 
efficients of the electric field excited by the incident light 
in the bulk part of the system. Ej Notice that Equation 
( p8| ) involves a selection rule on the spatial symmetry of 
the electric field. Since the transfer matrix Q commutes 
with the symmetry operations of the point group rele- 
vant to the incidence surface, the eigenstates of Q are 
classified according to the irreducible representations of 
that point group. Hence, c n must become zero if itbuik is 
attributed to a different irreducible representation from 
that of u\ . 

In order to confront our estimation with the predic- 
tions of group theory, let us consider again the previ- 
ous example of normal incidence along the (111) direc- 
tion. In this case incident light has a fcy = compo- 
nent within the surface Brillouin zone associated with 
the {111} set of planes. Following the abovementioncd 
procedure, we employed a 32-monolayer bare opal slab to 
obtain the electric field distribution in the space between 
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the 16th and 17th layers, where the bulk configuration 
is assumed to be achieved. At the same time, the pho- 
tonic dispersion relation of the infinite periodic system 
was also re-calculated for increasing values of k z . Figure 
5 shows both dispersion relation and symmetry assign- 
ment within the [1.10,1.3] interval. Now, we will pay par- 
ticular atcntion to the pair of frequency values marked 
with horizontal lines: for one of them (a/A = 1.15) no 
uncoupled bands are predicted by group theory, while 
for the other (a/ A = 1.25) some of that bands are ex- 
pected. In the first case, we obtained \c n \ 2 equal to 
0.689 and 0.348 for E bands with positive group veloc- 
ity and |c„| 2 = 0.007,0.006 for the ones with negative 
slope. When calculated for the higher frecuency, the val- 
ues of |c„| 2 for E(+) and E{-) bands are {0.265, 0.377, 
0.109, 0.035} and {0.085, 0.016, 0.098, 0.054} respec- 
tively. With regards to the uncoupled bands (i.e, those 
with A\ symmetry), their coefficients are less than 10 -16 . 
Therefore we can conclude that light transmission is ac- 
tually forbidden for such modes, as far as they do not 
contribute to the excited field inside the slab. This con- 
stitutes a confirmation of group theory predictions. 



VI. SUMMARY 

We have analyzed the symmetry properties of eigen- 
states along the high symmetry directions of close-packed 
bare opals according to the group theory. We found that 
some bands cannot be coupled with an external plane 
wave along several directions because of symmetry rea- 
sons, which was confirmed by layer-KKR calculations; 
nevertheless, the presence of uncoupled modes is not ex- 
pected to provide any observable features in the trans- 
mittance for this particular system. 
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satisfy Equation (|17|) but in an approximate way. How- 
ever, a suitable algorithm can always be found in order 
to estimate the character of computed eigenfunctions. 

In our present work, we have defined a strategy based 
on direct evaluation of Equation (^) , which will be now 
rewritten for the sake of clarity: 



P g = AE g = ^2 ' 



(Al) 



In the case of one-dimensional irreducible representations 
(g = 1), the character x(A) is given by 



Xid = a 



P E 

w 



(A2) 



Similar expressions can be obtained for two-dimensional 
(.9 < 2) representations 

(-lp+^Pj ■ E k )(E k ■ Ej) - (Pj ■ Ej)|E k | 2 ] 
X2D ^ IE 1 | 2 |E 2 | 2 -(E 1 -E 2 )(E 2 -E 1 ) 



(A3) 



and also for three-dimensional (g < 3) ones, given that 
some auxiliary quantities be defined: 



X3D 



W X • (v 2 X V 3 ) + V! • (w 2 X V 3 ) + Vx • (v 2 X w 3 ) 



vi ■ (v 2 x v 3 ) 



Ek ■ E; 



(v k )j = J^k ' 
(wk)j = P k Ej 



(A4) 

(A5) 
(A6) 



We finally have to mention that all symmetry assign- 
ments have been confirmed by means of a different nu- 
merical method for character evaluation which was pro- 
posed in a previous paper £9 
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APPENDIX A: NUMERICAL ASSIGNMENT OF 
IRREDUCIBLE REPRESENTATIONS 



APPENDIX B: COMPATIBILITY RELATIONS 
OF IRREDUCIBLE REPRESENTATIONS IN THE 
FCC LATTICE 

In Tables IV to IX, the compatibility relations between 
the irreducible representations for the points at the cor- 
ners of the irreducible part of the fee first Brillouin zone 
and those for the wave vectors on intermediate segments 
are presented. It can be observed that only one mirror 
reflection is defined on many of these segments because 
of their low symmetry invariance. 



As stated in Section III, the label assignment for a 
given mode E£ (r) is based on the evaluation of its char- 
acter under the symmetry operators within Gk- The first 
thing that needs to be remarked is that solutions ob- 
tained by means of the plane wave expansion method 
(in fact by means of any numerical calculation) do not 
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TABLE I. The symmetry point groups at the corner of the 
irreducible part of the Brillouin zone in the fee lattice accord- 
ing to Schoenflies notation. Notice that the components of 
representative vectors are in units of 27r/a 

Label Representative k vector Symmetry point group 



r 

X 
L 
U 

w 

K 



(0, 0, 0) 
(1, 0, 0) 
(1/2, 1/2, 1/2) 
(1, 1/4, 1/4) 

(1, 1/2, 0) 
(3/4, 3/4, 0) 



O h 
D 3d 

C2v 

D 2d 

C 2v 



28 A. Modinos, N. Stefanou and V. Yannopapas, Opt. Exp., 
8, 197 (2001) 

29 Further details can be found in T. Ochiai and J. Sanchez- 
Dehesa, Phys. Rev. B (in press) 
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Figure l.The first Brillouin zone of the fee lattice. 
The lines between white circles define the boundaries of 
the irreducible part of the zone. 

Figure 2. Photonic band structure along the XU and 
UL directions for a bare opal system. Labels refer to 
irreducible representations of C*2« and C\h respectively. 
An inset figure is included for the sake of clarity. 

Figure 3. Photonic band structure along the LT and 
TX directions for a bare opal system. Labels refer to 
irreducible representations of Cs v and C± v respectively. 
Some inset figures are included for the sake of clarity. 

Figure 4. Photonic band structure along the XW and 



WK directions for a bare opal system. Labels refer to 
irreducible representations of C2 V and C\h respectively. 
Some inset figures are included for the sake of clarity. 

Figure 5. Detail of the photonic band structure nor- 
mal to the {111} set of planes of a bare opal system. 
Notice that the maximum value of \k z \ is V3/2 in units 
of 27r/a. The thin horizontal lines define the frequencies 
at which the branching ratios are obtained. 
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TABLE II. The irreducible representations corresponding 
to eigenvalues first to twentieth at the corners of the irre- 
ducible part of the first Brillouin zone in a bare opal system 



r[Oh] 


X[D 4h ] 


L[D sd ] 


U[C 2v ] 


W[D 2d ] 


K[C 2v ] 


Tin 


E u 


E 9 


B 2 


Ai 


B 2 


Tig 


Eg 


E u 


Bi 


E 


Bi 


E u 


A 2g 


Ai„ 


Ai 


Bi 


Ai 


T 2u 


Eg 


E u 


A 2 


A 2 


A 2 


Eg 


Big 


E 9 


Bi 


E 


Bi 


Ti u 


Biu 


Alg 


B 2 


B 2 


B 2 


T 29 


E u 


E 9 


A 2 


A 2 


A 2 


T 2 u 


A 2u 


A 29 


B 2 


B 2 


B 2 




B 2 „ 


E u 


Ai 


E 


Ai 




Ai„ 


A 2u 


Bi 


Bi 


Bi 




Big 


A 29 


Bi 


Ai 


Bi 




Eu 


E 9 


B 2 


Bi 


B 2 




Eg 


E u 


A 2 


Ai 


A 2 




B 29 




A 2 


E 


A 2 








B 2 




B 2 








Bi 




Bi 








Ai 




Ai 








B 2 




B 2 








Ai 




Ai 








Bi 




Bi 



TABLE V. The compatibility relations of irreducible rep- 
resentations at X point in the fee lattice 



X[D 4h ] 


rx[c 4v ] 


XU[C 2v ] 


XW[C 2v ] 


Alg 


Ai 


Ai 


Ai 


A 2 g 


A 2 


B 2 


B 2 


Big 


Bi 


B 2 


Ai 


B 2 g 


B 2 


Ai 


B 2 


E 9 


E 


A 2 ffiBi 


A 2 0Bi 


Aiu 


A 2 


A 2 


A 2 


A 2 u 


Ax 


Bi 


Bi 


Biu 


B 2 


Bi 


A 2 


B 2 u 


Bi 


A 2 


Bi 


Eu 


E 


AiffiB 2 


Ai0B 2 



TABLE VI. The compatibility relations of irreducible rep- 
resentations at L point in the fee lattice 



L[D 3d ] 


TL[C 3v ] 


LW[C 2 ] 


LU[C lh ] 


LK[C lh ] 


Alg 


Ai 


A 


A' 


A' 


A 29 


A 2 


B 


A" 


A" 


E 9 


E 


A©B 


A'©A" 


A'©A" 


A lu 


A 2 


A 


A" 


A" 


A 2 u 


Ai 


B 


A' 


A' 


Eu 


E 


A©B 


A'eA" 


A'eA" 



TABLE III. The representations of incident plane waves 
and the subsequent uncoupled bands along some relevant 
high-symmetry directions in the Brillouin zone of the fee lat- 
tice. Point symmetry groups are described (in brackets) ac- 
cording to Schoenflies notation. 



Incidence Direction 


PW Representation 


Uncoupled Bands 


TL[C 3 u] 

rx[c 4v ] 

TK[C 2 v] 


E 
E 

Si ®B 2 


Ai,A 2 
A 1 ,A 2 ,B 1 ,B 2 
Ai,A 2 



TABLE IV. The compatibility relations of irreducible rep- 
resentations at r point in the fee lattice 





rx[c 4v ] 


ru[c lh ] 


TL[C 3v ] 


TK[C 2v ] 


rw[c lh ] 


Alg 


Ax 


A' 


Ai 


Ax 


A' 


A 29 


Bi 


A" 


A 2 


B 2 


A' 


Eg 


AiffiBi 


A'©A" 


E 


AieB 2 


2A' 


Tig 


E©A 2 


A'©2A" 


E©A 2 


A 2 eBieB 2 


A'©2A" 


T 2 g 


E©B 2 


2A'©A" 


EeAi 


AieA 2 eBi 


A'©2A" 


Aiu 


A 2 


A" 


A 2 


A 2 


A" 


A 2u 


B 2 


A' 


Ai 


Bi 


A" 


Eu 


A 2 ©B 2 


A'eA" 


E 


A 2 eBi 


2A" 


Tiu 


E©Ai 


2A'eA" 


EeAi 


AieBieB 2 


2A'eA" 


T 2 u 


E©Bi 


A'©2A" 


EeA 2 


AieA 2 eB 2 


2A'eA" 



TABLE VII. The compatibility relations of irreducible rep- 
resentations at U point in the fee lattice 



U[C 2v ] 


ru[c lh ] 


XU[C 2v ] 


LU[C lh ] 


UW[C lh ] 


Ai 


A' 


Ai 


A' 


A' 


A 2 


A" 


A 2 


A" 


A" 


Bi 


A' 


Bi 


A' 


A" 


B 2 


A" 


B 2 


A" 


A' 



TABLE VIII. The compatibility relations of irreducible 
representations at W point in the fee lattice 



W[D 2d ] 


rw[c lh ] 


XW[C 2v ] 


LW[C 2 ] 


UW[C lh ] 


WK[C lh ] 


Ai 


A' 


Ax 


A 


A' 


A' 


A 2 


A" 


A 2 


B 


A" 


A" 


Bi 


A" 


A 2 


A 


A" 


A" 


B 2 


A' 


Ax 


B 


A' 


A' 


E 


A'eA" 


Bx©B 2 


AeB 


A'eA" 


A'©A" 



TABLE IX. The compatibility relations of irreducible rep- 
resentations at K point in the fee lattice 



K[C' 2v ] 


TK[C 2v ] 


LK[C lh ] 


WK[C lh ] 


Ax 


Ax 


A' 


A' 


A 2 


A 2 


A" 


A" 


Bi 


Bi 


A' 


A" 


B 2 


B 2 


A" 


A' 
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